Environmental  correlates  of  Arctic  ice-edge  noise 

Nicholas  C.  Makrisa)  and  Ira  Dyer 

Department  of  Ocean  Engineering,  Massachusetts  Institute  of  Technology,  Cambridge, 

Massachusetts  02139 

( Received  5  February  1991;  accepted  for  publication  8  July  1 99 1 ) 

Temporal  variations  of  low-frequency,  broadband  ambient  noise  measured  in  early  summer 
under  drifting  ice  floes  of  the  marginal  ice  zone  (MIZ)  are  cross  correlated  with  local 
environmental  forces  and  ice  field  descriptors.  Surface  gravity  wave  forcing  is  the  primary 
correlate  of  the  noise;  its  interaction  with  ice  floes  generates  sound,  most  likely  via  flexural  floe 
failure,  and  unloading  motion,  within  a  few  kilometers  of  the  ice  edge.  Ice  concentration  is  also 
well  correlated  with  the  noise,  most  likely  parametrically.  That  is,  as  ice  concentration 
increases,  so  does  the  density  of  potential  sound  sources  in  the  ice  field.  Ice  stress,  ice  moment, 
and  wind  stress  magnitude,  while  highly  correlated  with  low-frequency  noise  in  the  fully  ice- 
covered  Arctic,  are  poor  correlates  in  the  MIZ.  Lateral  melt  rate,  as  a  surrogate  for  thermally 
induced  ice  stress,  is  also  poorly  correlated  with  low-frequency  MIZ  noise.  When  surface  wave 
forcing  is  weak,  about  1/2  the  time  in  this  experiment,  episodes  of  high  noise  are  sometimes 
observed.  These  episodes  are  roughly  3-6  h  in  duration  and  12  h  in  periodicity.  Evidence 
suggests  that  they  are  related  to  the  formation  and  advection  of  bands  of  highly  concentrated 
ice  by  internal  waves  during  off-ice  winds. 

PACS  numbers:  43.30.Ma,  43.30.Nb 


SYMBOLS  FOUND  IN  TABLES 

vg 

absolute  current  velocity  vector  at  50-m 

depth 

Symbols  used  specifically  in  the  tables  are  defined  here. 

Vr 

pressure  gradient  shear  stress  vector 

Subscripts  x,y  denote  east  and  north  components  respective- 

Mr 

lateral  melt  rate 

ly- 

SW 

surface  gravity  wave  amplitude 

^div 

ice  field  divergence  rate 

^12 

ice  field  shear  strain  rate 

1*1 

ambient  noise  rms  pressure  in  the  25-  to 

^xx  ^yy 

ice  field  east-west  expansion  rate 

50-Hz  band 

c 

ice  field  concentration 

Ta 

wind  stress  vector 

^ xx  *yy 

ice  field  east-west  expansion 

Tw 

current  stress  vector 

^abs 

ice  field  rotation  rate 

K 

ice  drift  velocity  vector 

(£div)2+(*u)2 

ice  field  deformation  energy 

INTRODUCTION 

As  part  of  the  marginal  ice  zone  experiment  of  1984 
(MIZEX  84),  MIT  measured  ambient  noise  beneath  freely 
drifting  ice  floes.  These  measurements  were  made  in  the 
Fram  Strait  near  80.5°  N  and  8°  E,  during  the  month  of  June. 
Our  purpose  in  acquiring  the  data  is  to  gain  an  understand¬ 
ing  of  the  physical  mechanisms  responsible  for  marginal  ice 
zone  (MIZ)  underwater  noise. 

In  this  paper,  we  compare  temporal  variations  of  low- 
frequency  ambient  noise  with  environmental  forces  applied 
to  the  ice,  including  those  due  to  wind,  current,  drift,  surface 
gravity  wave  amplitude  and  heat  flux.  We  find  that  over  long 
periods,  about  3  h  or  more  (autocorrelation  e-folding  time), 
low-frequency  MIZ  noise  cross  correlates  best  with  surface 
gravity  wave  amplitude.  Surface  gravity  waves  penetrate 
into  the  MIZ  from  the  open  waters  of  the  East  Greenland 
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Sea,  and  are  primarily  responsible  for  the  break-up  of  ice 
floes. 

We  also  compare  temporal  variations  of  the  noise  with 
certain  ice  descriptors.  These  include  ice  concentration,  de¬ 
formation,  and  deformation  rate.  We  find  that,  over  long 
periods,  the  noise  is  well  correlated  with  ice  concentration. 
For  the  shorter  periods,  and  when  surface  wave  forcing  is 
weak,  episodes  of  high  noise  are  sometimes  observed.  These 
episodes  are  of  1-  to  2-h  autocorrelation  time  ( 3-  to  6-h  dura¬ 
tion)  and  12-h  periodicity.  Qualitative  evidence  suggests 
that  they  are  related  to  the  formation  and  advection  of  bands 
of  highly  concentrated  ice  by  internal  waves.  Together,  these 
findings  suggest  that  low-frequency  ambient  noise  is  para¬ 
metrically  dependent  upon  ice  concentration  in  the  MIZ, 
over  both  short  and  long  periods.  In  particular,  we  suggest 
that  as  ice  concentration  increases,  so  does  the  density  of 
potential  surface  sound  sources. 

Over  long  periods  internal  ice  stress,  applied  moment, 
and  the  individual  stresses  resulting  from  wind,  current,  and 
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drift  are  generally  poor  correlates  of  noise  in  the  MIZ.  This  is 
in  contrast  to  the  results  of  Makris  and  Dyer1  for  pack-ice- 
covered  waters  of  the  central  Arctic.  We  interpret  this  differ¬ 
ence  to  be  a  direct  consequence  of  the  reduced  fetch  available 
for  normal  horizontal  stress  and  moment  accumulation  in 
the  MIZ.  That  is,  decreased  ice  floe  size  and  less  than  contin¬ 
uous  ice  cover,  as  well  as  the  proximity  of  the  ice  edge,  where 
normal  horizontal  stresses  may  be  relieved,  are  primarily 
responsible  for  the  reduction  of  ice  stress  in  the  MIZ  with 
respect  to  the  central  Arctic. 

Diachok  and  Winokur2  first  showed  that  ambient  noise 
of  the  MIZ  could  vary  from  that  of  the  adjacent  open  ocean 
and  pack-ice-covered  seas.  Via  a  series  of  across  ice-edge 
synoptic  (2  h)  sonobuoy  measurements,  they  observed  in¬ 
creased  noise  levels  at  or  near  the  ice  edge,  which  they  attri¬ 
buted  to  surface  gravity  wave  interaction  with  ice  floes. 
From  measurements  made  via  a  towed  array,  within  a  1-day 
period,  Yang  et  al.3  more  recently  found  the  along  ice-edge 
distribution  of  noise  to  be  localized  at  “hotspots”  of  less  than 
5-km  length  and  about  50-km  spacing.  They  speculated  that 
ice-edge  eddies  may  have  been  responsible  for  the  hotspots, 
although  they  lacked  environmental  data  to  support  this. 
Subsequently,  Johannessen  et  al.4  synoptically  measured 
ambient  noise  levels  with  sonobuoys,  within  a  6-h  period,  in 
a  region  containing  two  mesoscale  eddies,  and  found  indica¬ 
tions  of  hotspots,  although  their  spatial  sampling  interval 
was  at  roughly  five  times  the  reported  hotspot  size.  In  con¬ 
trast,  the  present  study  has  longer  noise  records,  finer  spatial 
sampling  of  the  noise  field  and  a  wider  range  of  environmen¬ 
tal  data,  but  it  was  carried  out  in  a  region  with  no  apparent 
eddy  motion,  and  did  not  include  observations  either  in  the 
adjacent  open  water  or  fully  covered  central  Arctic  ice  field. 

Lewis  and  Denner5  have  studied  the  relationship 
between  noise  and  heat  flux.  They  have  found  that  episodes 
of  high  noise  in  the  frequency  range  up  to  several  kHz  corre¬ 
late  well  with  heat  flux.  This  mechanism  is  thought  to  be 
important  when  thermal  changes  are  large,  as  is  the  case  in 
the  Arctic  for  most  of  the  year,  but  typically  not  in  the  sum¬ 
mer.  Because  of  the  summer  season,  and  the  lower  frequency 
range  considered,  the  MIZEX  84  noise  data  we  have  studied 
are  unlikely  to  relate  to  heat  flux.  But  it  is,  nonetheless,  in¬ 
cluded  here  as  a  possible  environmental  correlate. 

Acoustic  measurements  were  made  using  hydrophones 
suspended  at  60  m  below  freely  drifting  ice  floes.  While  up  to 


24  such  hydrophones  were  deployed,  roughly  half  this 
number  were  within  our  operational  rf  telemetry  range  at  a 
given  time.  The  operational  phones  formed  a  typically  5-km 
scale  random  array.  In  response  to  a  wide  range  of  wind  and 
current  magnitudes  and  directions,  this  array  drifted  within 
5-20  km  of  relatively  open  water.  Via  an  integral  sensor 
tracking  system  (STS),6,7  floes  equipped  with  acoustic  sen¬ 
sors  and  current  meters  could  be  located.  These  locations 
were  also  used  to  calculate  ice  field  descriptors  such  as  con¬ 
centration  and  strain.  Wind  and  drift  measurements  were 
made  from  the  experiment  base  ship  Kvitbjom  (Kv),  which 
was  part  of  the  STS.  While  the  ice-edge  position,  roughly 
normal  to  the  concentration  gradient,  varied  by  tens  of  kilo¬ 
meters  over  the  course  of  the  roughly  2  week  acoustic  experi¬ 
ment,  it  maintained  an  east-west  orientation  within  20  km  of 
the  Kv.  Charts  of  geostrophic  current  and  microwave  ima¬ 
gery  show  that  there  were  no  prominent  mesoscale  eddies  in 
this  vicinity  at  any  time  during  the  experiment.8-10  This  was 
confirmed  via  ice  kinematic  data.11 

Noise  in  a  25-  to  50-Hz  band  is  used  as  a  surrogate  for 
the  broad  distribution  of  low-frequency  ambient  noise,  as 
may  be  justified  by  inspection  of  Fig.  1,  which  indicates  that 
this  band  is  at  or  near  the  peak  of  the  noise  spectrum.  Spectra 
from  about  5  to  300  Hz  at  other  times  increased  or  decreased 
nearly  uniformly,  so  that  the  octave  containing  the  peak  clo¬ 
sely  mimics  temporal  variations  of  the  low-frequency  spec¬ 
trum. 

Records  of  rms  noise  in  the  25-  to  50-Hz  band  were 
constructed  from  the  analog  and  digital  recordings  of  each 
available  drifting  sensor.  The  noise  records  of  two  of  the 
most  consistently  operational  sensors  appear  in  Fig.  2,  and 
were  obtained  by  integrating  and  then  square  rooting  the 
pressure  spectral  density  over  the  25-  to  50-Hz  band.  The 
averaging  period  for  each  sample  varied  from  1  to  5  min.  The 
sampling  interval  of  these  long-term  noise  records  was  no 
greater  than  1/2  h,  and  sometimes  finer.  Gaps  in  the  ambient 
noise  records  were  caused  by  contamination  from  interfering 
acoustic  experiments,  and  intermissions  in  recording,  the 
statistics  of  which  appear  in  Table  I.  Due  to  a  long,  roughly 
3-day  intermission  in  recording,  the  noise  records  occurring 
before  and  after  this  intermission  have  been  analyzed  separa¬ 
tely.  They  are  respectively  referred  to  as  the  first  and  second 
drift  phases. 

By  comparing  the  noise  records  of  the  two  sensors 


FIG.  I.  Composite  of  ambient  noise  ob¬ 
served  on  Julian  Day  177.4  with  sensor 
D2.  Data  were  taken  sequentially  over 
several  decades  in  frequency,  and  then 
repeated  to  cover  adjacent  decadal 
ranges.  All  data  were  acquired  within 
1/2  h.  Peaks  at  about  3  and  6  Hz  are 
caused  by  hydrophone  cable  strum.  The 
broad  peak  from  about  5  to  7  kHz  is  an 
artifact  associated  with  the  cable/ 
preamplifier  electronics,  and  should  be 
ignored.  A  40-dB/decade  high-pass 
filter  with  comer  frequency  at/=  1  Hz 
has  whitened  the  noise  for/<  1  Hz,  and 
should  be  corrected  for  upon  use  of  this 
graph  in  this  frequency  range. 
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FIG.  2.  Root-mean-square  noise  in  the  25-  to  50-Hz  band  for  sensors  ( a )  D6 
and  (b)  D24.  These  sensors  were  separated  by  between  1.5  to  5.0  km 
throughout  the  time  period  shown.  Pluses  indicate  data  points  within  each 
of  two  separate  drift  phases.  The  first  drift  phase,  from  J 168  to  J 174  has  a 
mean  of  0.07  and  0.06  Pa  and  a  standard  deviation  of  0.03  and  0.02  Pa  for 
(a)  and  (b),  respectively.  The  second  drift  phase,  from  J 177  to  J 181,  has  a 
mean  of  0. 1 5  and  0. 1 2  Pa  and  a  standard  deviation  of  0. 1 2  and  0.09  Pa  for 
(a)  and  (b),  respectively. 

shown  in  Fig.  2,  it  is  evident  that  long-term  temporal  varia¬ 
tions  of  the  noise  are  nearly  uniform  over  the  roughly  5-km 
region  occupied  by  the  acoustic  array.  This  observation  is 
quantitatively  confirmed  via  cross  correlation  of  the  respec¬ 
tive  noise  records.  A  peak  correlation  coefficient  of  0.79  was 
obtained  during  the  first  drift  phase  and  0.97  during  the  se¬ 
cond,  each  at  zero  time  lag.  The  cross-correlation  algorithm 
that  we  used  required  no  interpolation  of  data  gaps  and  so 
introduced  no  hypothetical  information.  This  algorithm,  de¬ 
scribed  in  the  Appendix,  was  used  for  all  auto  and  cross 


TABLE  I.  Gap  statistics  for  long-term  noise  records. 


Record 

sensor  drift  phase 

Mean  gap 
(h) 

Longest  gap 
(h) 

%  gap 

D6 

1st 

2.3 

10.4 

41 

D24 

1st 

3.3 

14.4 

53 

D6 

2nd 

2.3 

7.0 

34 

D24 

2nd 

2.3 

7.0 

34 

correlations  in  this  paper.  Because  the  noise  record  for  sen¬ 
sor  D6  was  better  sampled  during  the  first  drift  phase,  it  is 
used  exclusively  in  this  paper  for  correlations  with  environ¬ 
mental  records. 

More  densely  sampled  acoustic  data  measured  during  a 
particular  high  noise  episode,  from  an  array  consisting  of 
several  hydrophones,  are  analyzed  in  Sec.  III. 


I.  NOISE  AUTOCORRELATION  TIME  IS  SHORTER  IN 
THE  MIZ 

Time  series  of  the  rms  noise,  as  shown  in  Fig.  2(a)  have 
been  autocorrelated  for  each  of  the  drift  phases,  and  are 
shown  in  Fig.  3.  Spectra  are  also  given.  The  e-folding  time 
(the  time  at  which  the  normalized  autocorrelation  falls  to 
e  ~ 1 )  for  the  two  are  listed  in  Table  II,  at  3  and  1 3  h,  respecti¬ 
vely.  These  are  considerably  shorter  than  observed  for  the 
central  Arctic,  about  29  h.1  Therefore,  the  available  data 
indicate  that  the  MIZ  is  an  acoustically  more  dynamic  re¬ 
gion  which,  of  course,  is  a  consequence  of  the  greater  dyna¬ 
mism  of  the  MIZ  environment. 

Spectral  slopes  shown  in  Fig.  3  for  large  values  of  v  are 
about  v-2,  and  compare  with  v~ 3/2  for  the  central  Arctic.1 
No  strong  semidiurnal  or  diurnal  components  are  observed 
in  either,  although  an  argument  can  be  made  for  a  weak 
semidiurnal  component  during  the  first  drift  phase.  The  fas¬ 
ter  evolution  of  the  first  drift  phase,  compared  to  the  second, 
shows  in  its  spectrum  as  well  as  in  its  autocorrelation  time. 

While  the  spectra  have  no  strong  semidiurnal  compon¬ 
ents,  we  will  show  in  Sec.  Ill  that  such  components  can 
sometimes  be  observed.  What  is  required  is  the  absence  of 


FIG.  3.  The  time  lagged  normalized  autocorre¬ 
lation  and  spectral  level  of  the  25-  to  50-Hz  rms 
noise  as  measured  by  sensor  D6  for  (a)  the  first 
drift  phase,  (b)  the  second  drift  phase.  For  the 
first  drift  phase,  a  3-h  e-folding  time  is  measured 
and  a  spectral  peak  at  v  =  2  cycle/day  can  be 
argued  for  in  the  spectrum,  corresponding  to 
the  oceanographically  important  period  of  iner¬ 
tial  oscillations.  For  the  second  drift  phase  a  13- 
h  e-folding  time  is  measured  and  there  is  no  evi¬ 
dence  of  a  spectral  peak  at  2  cycle/day. 
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surface  wave  forcing,  the  major  factor  in  noise  generation 
(Sec.  II),  and  other  environmental  effects  that  concentrate 
the  ice  (Sec.  III). 

II.  SURFACE  GRAVITY  WAVES  CORRELATE  BEST 
WITH  THE  NOISE 

A.  Environmental  forcing  functions 

A  variety  of  dynamical  forcing  functions  were  cross  cor¬ 
related  with  the  noise.  These  include  horizontal  wind  shear 
stress  which  acts  on  the  upper  surface  of  the  ice,  current 
shear  stress  which  acts  on  the  lower  surface,  normal  stress 
applied  by  the  Coriolis,  and  lineal  acceleration  forces  which 
act  uniformly  on  an  ice  floe’s  submerged  section,  lateral  melt 
rate  as  a  surrogate  for  heat  flux,  and  surface  gravity  wave 
amplitude.  Each  of  these,  individually  or  in  combination, 
can  cause  ice  fracture  and,  hence,  noise. 

The  following  expressions  describe  the  horizontal 
stresses  acting  upon  a  drifting  ice  floe,  as  induced  by  wind, 
current,  and  drift.  These  are  similar  to  those  found  in  Makris 
and  Dyer.1  In  particular,  wind  speed  |Va|  measured  on 
board  the  Kv  10  m  above  the  ice  was  averaged  over  and 
sampled  at  10-min  intervals,  and  converted  to  wind  shear 
stress  via 

Ta  ^10 Pa  I  y a  I  ^  a  > 

where  pa  is  air  density  and  cl0  =  0.0023  is  the  air-ice  drag 
coefficient  corresponding  to  neutrally  stable  conditions  en¬ 
countered  during  MIZEX  84. 12  Current  velocity  relative  to 
the  drifting  ice  measured  at  10-m  depth  and  averaged 
over  and  sampled  at  10-min  intervals,  were  converted  to  cur¬ 
rent  stress  magnitude  via 

K|  =cwpw\\ j2, 

where  cw  =  0.0046  is  the  ice-ocean  drag  coefficient  based 
upon  measurements  at  10-m  depth,  and  pw  is  the  density  of 
sea  water.  The  ice-ocean  shear  stress  is  rotated  clockwise 
downward  through  the  oceanic  boundary  layer  with  inflow 
angle 

cos  a  —  r  w'Vw/\tw  1 1 |,  a  =  -28°, 

where  the  numerical  value  has  been  determined  by  Mak¬ 
ris.11 

Coriolis  normal  stress  is  derived  from 

<rc  =  -2/olAn(sin^)(iXVf), 

where  V,  is  the  ice  drift  velocity,  pt  is  the  ice  density,  h  is  the 
ice  sheet  thickness,  taken  to  be  1.5  m,  ^  is  the  latitude,  fl  is 
the  Earth’s  angular  speed,  and  i  is  a  unit  vector  in  the  verti¬ 
cal.  Ice  drift  velocity  was  obtained  by  Kalman  filtering  satel¬ 
lite  navigation  data  for  the  Kv  with  a  6-h  low-pass  filter  to 
hourly  values13  which  were  later  interpolated  to  half  hourly 
values  by  us.  The  normal  stress  associated  with  lineal  acce¬ 
leration,  pih(d\i/dt),  was  determined  by  differentiating 
the  low-pass  drift  velocity.  The  normal  stress  applied  by  the 
pressure  gradient  force  is 

d\„ 

crp  =  2plt)/isn(sin  <f>)  (iXVg)  +  pwhs  — , 

dt 

where  hs  is  the  submerged  thickness  of  the  ice  sheet  and  \g  is 
the  absolute  current  velocity.  It  is  measured  at  50-m  depth 


where  momentum  transfer  from  the  Ekman  layer  is  assumed 
to  be  negligible.  The  lineal  acceleration  term  is  of  the  same 
order  as  the  geostrophic  term  due  to  rapidly  varying  currents 
in  the  MIZ. 

The  internal  ice  stress  acting  normal  to  the  ice  sheet’s 
vertical  section  is  approximated  as1 

S=  (ra  +t  1,)(Ls/h) 

where  the  combined  contributions  of  Coriolis,  pressure  gra¬ 
dient,  and  lineal  acceleration  normal  horizontal  stresses  are 
negligible  for  our  experiment.  The  stress  moment  acting 
about  the  ice  sheet’s  central  horizontal  plane  is  approximat¬ 
ed  as1 

M  =  ix  (ra  +  rw )  (Lm/ 2). 

Lengths  parallel  to  S  and  normal  to  M  are  Ls  and  Lm , 
respectively,  through  which  horizontal  load  and  bending 
moment  may  be  accumulated.  We  estimate  Ls  to  be  less  than 
or  equal  to  a  few  floe  lengths  or  no  more  than  about  500  m, 
and  Lm  to  be  less  than  or  equal  to  a  typical  floe  length  or  no 
more  than  about  100  m. 

Lateral  melt  rate,  Mr,  is  assumed  to  be  roughly  propor¬ 
tional  to  AT i;36, 14  where  A Tw  =  Tw  —  Tf  is  the  elevation  of 
sea  temperature  above  freezing  at  Tf—  —  1.98  °C.  Sea  tem¬ 
perature  Tw  was  measured  at  5  m  below  sea  level,  and  aver¬ 
aged  and  sampled  at  10-min  intervals. 13  For  our  experiment, 
Mr  is  a  surrogate  for  heat  flux. 

We  computed  the  surface  gravity  wave  field  in  the  East 
Greenland  Sea  adjacent  to  the  ice  edge  from  contemporary 
meteorological  charts  via  the  Pierson  et  al  spectral  meth¬ 
od.  15  This  proved  to  be  quite  simple  since  only  on-ice  winds 
have  sufficient  fetches  to  generate  potentially  important  sur¬ 
face  wave  amplitudes  at  the  ice  edge.  Surface  wave  ampli¬ 
tudes  within  the  MIZ,  at  the  acoustic  array,  were  determined 
for  wave  periods  in  excess  of  8  s  from  empirically  and  theo¬ 
retically  derived  attenuation  coefficients,16  and  estimates  of 
the  array’s  distance  from  the  ice  edge.  Waves  in  this  period 
regime  are  subsequently  referred  to  as  swell,  at  least  because 
wind  stress  in  the  ice  field  no  longer  drives  the  gravity  waves. 
We  did  not  consider  waves  with  periods  shorter  than  8  s 
because  investigations  of  flexural  floe  failure  by  Goodman  et 
al.;11  and  high-frequency  floe  collisions  ( ~  1  Hz)  in  the  MIZ 
by  Martin  and  Becker18  indicate  that  surface  waves  in  the 
swell  regime  (periods  >8  s)  are  primarily  responsible.  In 
addition,  earlier  work  by  Diachok  indicates  the  higher  the 
wave  period,  the  more  intense  the  low-frequency  noise.19 

The  e-folding  autocorrelation  times  of  the  ice-forcing 
functions  are  summarized  in  Table  II.  During  the  first  drift 
phase,  the  roughly  3-h  time  scale  of  the  noise  was  most  simi¬ 
lar  to  time  scales  of  stresses  derived  from  ice  drift  and  cur¬ 
rent.  Of  these,  only  drift  speed  (Fig.  4)  and  current  stress 
(not  shown)  display  evidence  of  semidiurnal  peaks  in  their 
spectra,  as  does  the  noise,  the  latter  only  weakly.  During  the 
second  drift  phase,  the  e-folding  time  of  the  noise  increased 
by  nearly  a  factor  of  4.  It  is  much  shorter  than  the  e-folding 
time  calculated  for  surface  waves,  but  longer  than  that  mea¬ 
sured  for  wind  stress  magnitude. 

Results  of  cross  correlation  between  dynamical  forcing 
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TABLE  II.  Autocorrelation  e-folding  times  for  noise  and  environmental 
forcing  functions,  in  hours.  (Data  sampled  at  1/2-h  intervals. )  Incomplete 
environmental  data  in  the  second  drift  phase  account  for  the  lack  of  some 
entries. 


1st  drift 

2nd  drift 

\p\ 

3 

13 

kal 

8 

6 

kJ 

2 

K  +  rj 

5 

ka  -Tw\ 

4 

|V,| 

4 

4 

\d\x/dt  | 

3 

3 

l°>l 

3 

Mr 

19 

SW 

23 

48 

functions  and  noise  are  summarized  in  Table  III.  During 
both  drift  phases,  surface  gravity  wave  amplitude  is  by  far 
the  best  correlate  of  the  noise,  supporting  this  as  a  general 
result.  The  lower  correlation  found  during  the  first  drift 
phase  is  primarily  due  to  episodes  of  high  noise  which  occur 
during  weak  surface  wave  activity.  Evidently,  these  noise 
episodes  are  not  caused  by  surface  waves.  They  will  be  dis¬ 
cussed  in  Sec.  III. 

Figures  4  and  5,  for  the  first  and  second  drift  phases 
respectively,  show  time  series  for  measured  noise  and  calcu¬ 
lated  surface  wave  amplitude.  Visual  comparison  of  the  two 
supports  the  cross-correlation  results  of  Table  III,  from 
which  we  have  concluded  that  surface  wave  correlation  is 
the  best.  Another  environmental  force  that  has  scored  well 
in  Table  III  (wind  stress)  is  also  shown,  as  is  ice  concentra¬ 


TABLE  III.  Cross  correlations  between  ambient  noise  (rms  pressure)  in 
the  25-  to  50-Hz  band  and  environmental  forcing  functions.  Environmental 
time  series  are  held  fixed.  (Data  sampled  at  1/2-h  intervals. )  [Peak  correla¬ 
tion  coefficient,  time  lag  (h).] 


1st  drift 

\P\ 

2nd  drift 

\P\ 

K\ 

0.15,  -  1.5 

0.56,  0 

kJ 

0.35,  -0.5 

ka  +  Tw\ 

0.25,  -  1.0 

\Ta-rw\ 

0.23,  -0.5 

|V,| 

0.13,  1.0 

-0.27,  2.6 

\d\{/dt\ 

WP\ 

0.23,  -2.5 

Mr 

0.26,  3.0 

SW 

0.93,  0 

tion  (to  be  discussed  subsequently).  Arguments  in  favor  of 
these  clearly  can  be  made,  and  we  will  do  so  for  ice  concen¬ 
tration,  but  wind  stress  is  a  far  less  consistent  correlate,  and 
does  come  out  second  best  to  surface  waves. 

The  primary  cause  for  the  break-up  of  ice  floes  in  the 
MIZ  is  surface-wave-induced  flexural  failure,  as  may  be  sur¬ 
mised  from  the  work  of  Goodman  et  at.17  Their  results  are 
the  basis  for  Fig.  6,  which  gives  the  surface  wave  amplitude 
required  for  flexural  failure.  Comparing  Figs.  4  and  5  with 
Fig.  6,  we  observe  that  surface  gravity  wave  amplitudes  were 
sufficient  to  cause  flexural  floe  failure  within  several  kilo¬ 
meters  of  the  ice  edge  during  at  least  some  periods  of  high 
noise  within  each  drift  phase.  Low-frequency  sound  then 
could  have  been  radiated  into  the  water  via  the  fracture  itself 
and  from  unloading  motion  of  ice  floes  once  they  failed.20 


FIG.  4.  Noise  and  various  environmental 
records  during  the  first  drift  phase:  (a) 
Noise  as  registered  by  D6.  (b)  Ice  concen¬ 
tration  derived  from  the  STS  strain  array. 
Initializations,  scaling  and  interpolations 
were  based  upon  ice  concentration  obtained 
from  remote  images,  (c)  Significant  wave 
amplitude,  S  W,  for  surface  gravity  waves  in¬ 
cident  upon  the  ice  edge  with  period  >8  s. 
Significant  wave  amplitude  is  the  amplitude 
of  the  1/3  highest  waves. 15  To  obtain  SW  in 
the  MIZ  at  the  acoustic  array,  at  a  maxi¬ 
mum  penetration  of  20  km  into  the  ice,  mul¬ 
tiply  the  left-hand  scale  by  0.05.  To  obtain 
SW  for  the  full  incident  spectrum,  multiply 
the  left-hand  scale  by  18.  Surface  wave  am¬ 
plitude  relates  to  wind  stress  and  wind  direc¬ 
tion  with  a  time  lag  the  order  of  1  day.  (d) 
Wind  stress  magnitude;  see  text  for  conver¬ 
sion  from  wind  speed,  (e)  Wind  vector  dir¬ 
ection,  clockwise  from  true  north.  The  dir¬ 
ection  given  is  the  one  the  wind  is  pointed 
toward,  (f)  Ice  drift  speed  magnitude. 
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FIG.  5.  Noise  and  various  environmental  re¬ 
cords  during  the  second  drift  phase:  (a)  Noise 
as  registered  by  D6.  (b)  Ice  concentration  der¬ 
ived  from  remote  images  and  field  observations, 
(c)  Significant  wave  amplitude,  SW,  for  surface 
gravity  waves  incident  upon  the  ice-edge  with 
period  >8  s.  To  obtain  SW  in  the  MIZ  at  the 
acoustic  array,  at  a  maximum  penetration  of  5 
km  into  the  ice,  multiply  the  left-hand  scale  by 
0.5.  To  obtain  SW  for  the  full  incident  spec¬ 
trum,  multiply  the  left-hand  scale  by  8.  Surface 
wave  amplitude  relates  to  wind  stress  and  wind 
direction  with  a  time  lag  the  order  of  1  day.  (d) 
Wind  stress  magnitude,  (e)  Wind  vector  direc¬ 
tion,  clockwise  from  true  north.  The  direction 
given  is  the  one  the  wind  is  pointed  toward.  (0 
Ice  drift  speed  magnitude. 


The  center  frequency  of  fracture  radiation  is  roughly 
/  =  cr/L ,  where  cr  is  the  Rayleigh  wave  speed  ( s  1700  m/s) 
which,  for  long  fracture  lengths  L  ( >35  m)  is  in  the  octave 
band  considered.  Also,  unloading  motion  radiates  at 
/  =  cg/xe,  where  cg  is  the  group  speed  of  flexural  waves  in 
the  ice  floe  and  xe  is  the  distortion  distance;  this  lies  within 
the  range  pertinent  to  our  study,  10  </ <  100  Hz. 


FLOE  SIZE  ,  m 


FIG.  6.  Minimum  significant  wave  amplitude  needed  to  cause  flexural  fai¬ 
lure  in  a  floe  for  the  given  ocean  surface  wavelength  A.  and  floe  length.  The 
values  are  based  upon  a  model  by  Goodman  et  a/.17  where  the  wave  re¬ 
sponds  to  an  initially  rigid  floe,  upon  which  the  subsequent  hydrodynamic 
pressure  field  deforms  the  ice  elastically.  The  ice  thickness  is  taken  to  be  1.5 
m,  a  critical  strain,  ecrjt  =  4.3 X 10 ~ 5  and  a  pre-existing  crack  length  of  10 
mm  are  assumed  for  failure,  based  upon  observations  of  real  floes  in  swell. 1 7 
To  break  MIZ  floes,  typically  of  roughly  100-m  length,  the  least  surface 
wave  amplitude  required  is  in  the  swell  regime,  with  wavelengths  of  roughly 
100  to  150  m. 
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However,  it  is  also  possible  that  low-frequency  sound  is 
generated  by  ice  floes  in  relative  motion,  as  their  edges  tap 
against  each  other.  Both  horizontal  and  vertical  high-fre¬ 
quency  floe  collisions  in  the  MIZ  are  correlated  with  local 
swell  for  amplitudes  in  excess  of  <5—0.001  m,18  the  brine 
layer  spacing  of  sea  ice.21  For  smaller  wave  amplitudes, 
these  edge  asperities  apparently  cause  adjacent  floes  to  stick 
together.  We  think  of  8  as  the  correlation  length,  in  the  slide 
plane,  of  the  asperity  roughness.  Vibration  could  then  be 
caused  as  these  asperities  slide  past  each  other  at  frequency 
/—SW<ysw  /<5,  and  with  <ysw  <  1  rad/s,  it  can  be  seen  that /is 
in  the  low-frequency  range  of  interest. 

Finally,  it  might  be  thought  possible  that  sound  is  gener¬ 
ated  for  those  larger  relative  edge  motions  that  cause  the 
asperities  to  break.  The  frequency  estimate  for  this  is  ca/l, 
where  ca  is  the  group  speed  of  transverse  motion  of  an  asper¬ 
ity  of  roughness  height  /  and  thickness  8.  For  any  reasonable 
ratio  of  1/8,  the  frequency  range  of  such  radiation  is  much 
too  high  to  affect  our  observations. 

Considerable  evidence  to  date  strongly  supports  the  use 
of  center  frequency  as  a  diagnostic  for  the  source  mechan¬ 
ism.  The  reason  simply  is  that  an  individual  fracture,  colli¬ 
sion,  asperity  collapse,  etc.,  is  known  to  occupy  a  bandwidth 
about  equal  to  the  center  frequency,20,22  i.e.,  about  an  octave 
in  width.  As  a  relatively  narrow-band  transient,  we  are  thus 
justified  in  using  its  estimated  center  frequency  to  help  select 
the  more  plausible  mechanisms  from  all  those  that  can  be 
logically  proposed  for  the  observation  band  considered  in 
this  paper. 

B.  Ice  field  descriptors 

Ambient  noise  was  also  cross  correlated  with  a  variety 
of  ice  field  descriptors.  These  include  ice  concentration,  ex¬ 
pansion,  expansion  rate,  divergence  rate,  principal  shear 
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strain  rate,  rotation  rate,  and  deformation  energy.  These  de¬ 
scriptors  were  derived  from  the  strain  rate  tensor  of  the  ice 
field  which  was  calculated  via  a  linear  least-squares 
method23  from  the  relative  positions  and  velocities  of  eight 
freely  drifting  sensors  in  the  STS  network.  (Position  data  in 
the  STS  network  were  obtained  from  raw  intersensor  time 
delays,  hand  edited,  linear  filtered,  and  finally  averaged  over 
and  sampled  at  6-min  intervals. )  This  method  finds  the  velo¬ 
city  gradient  which  minimizes  the  mean  square  error 
between  measured  velocities  of  the  sensor  cluster,  and  those 
possible  in  a  linear  velocity  field.  The  strain  rate  tensor  is 
defined  as 

e..=±|^  +  ^A 

u  2  \dxj  dxj  ’ 

where  i  =  a,b  and  j  =  a,b\  and  where  xa,  xb  and  ua,ub,  re¬ 
spectively,  are  orthogonal  positions  and  velocities  in  a  trans¬ 
lating,  rotating,  horizontal  Cartesian  frame.  The  error  in  the 
strain  rate  tensor  due  to  the  linearization  of  the  velocity  field 
is  no  more  than  30%.  The  vorticity  of  the  cluster  determined 
within  the  Cartesian  frame  is 


(  dub 

dua  \ 

\dxa  " 

dxb  ) 

The  divergence  rate 

^div  ^aa  ^ bb  » 

principal  shear  strain  rate 


and  deformation  energy 

(*div)2+  (^12  )2 

are  invariant  to  rotations  of  coordinates.  The  east-west  ex¬ 
pansion  rate  (i.e.,  north-south  across  ice  contraction  rate) 

^ xx  ^yy* 

where  x  and  y  are  east  and  north  position  coordinates,  de¬ 
scribes  the  along-ice-edge  expansion  rate  of  the  ice  field  since 
the  orientation  of  the  ice  edge  is  essentially  east-west. 

The  temporal  variations  of  ice  concentration  and  east- 
west  expansion  were  determined  by  time  integration  of  the 
negative  divergence  rate  and  east-west  expansion  rate,  re¬ 
spectively.  These  estimates  are  presumably  stable  since  these 
descriptors  were  obtained  from  differentiation  of  the  raw 
position  data  and  are  consistent  with  estimates  using  subsets 
of  the  original  STS  array.  The  absolute  scale  of  the  ice  con¬ 
centration  estimate  as  well  as  interpolations  made  between 
J171.5  to  J174  and  J177  to  J181  were  determined  from  re¬ 
mote  images  of  the  ice  edge, 10,24,25  and  field  observations 
made  by  MIT. 

The  e-folding  autocorrelation  times  of  the  ice  field  de¬ 
scriptors  are  shown  in  Table  IV.  Divergence  rate,  shear 
strain  rate,  expansion  rate,  and  vorticity  all  have  e-folding 
times  substantially  shorter  than  that  of  the  noise.  These  time 
scales  reflect  the  heterogeneous  adjustments  of  individual 
floes  in  the  ice  field  as  they  respond  to  environmental  forc¬ 
ing.  Ice  concentration  and  expansion  rate  have  longer  time 
scales,  of  the  same  order  as  that  of  wind  direction.  The  longer 
time  scales  reflect  the  homogeneous  adjustments  of  the  ice 


TABLE  IV.  Autocorrelation  e-folding  times  for  ice  field  descriptors,  in 
hours.  (Data  sampled  at  6-min  intervals.) 


^div 

^12 

^xx  ^yy 

c 


^xx  ^yy 
*>abs 

KJ 

(^div  )2  +  (^12  )2 


1st  drift  2nd  drift 


1 

2 

1 

17  34 

89 
1 
2 
1 


field  to  variations  in  wind  direction  and  magnitude. 

The  noise  records  were  interpolated  to  a  6-min  sampling 
period  to  match  the  higher  sampling  of  the  ice  field  descrip¬ 
tor  records.  Results  of  cross  correlation  between  ice  field 
descriptors  and  noise  are  summarized  in  Table  V.  During 
the  first  drift  phase,  ice  concentration  correlates  best  with 
the  noise,  but  the  degree  of  correlation  is  not  particularly 
high.  An  explanation  for  this  will  be  given  in  Sec.  Ill  since  it 
pertains  to  noise  episodes.  Strain  rate  parameters  are  all  poor 
correlates  of  the  noise.  This  is  expected,  since  if  the  noise  is 
related  to  floe  collisions,  strain  rates  would  be  high  during 
low  ice  concentration,  when  few  collisions  are  possible. 

High  correlation  between  ice  concentration  and  noise 
exists  during  the  second  drift  phase  (although  no  other  ice 
field  descriptors  are  available  due  to  shut-down  of  the  STS 
network).  The  consistency  of  best  correlation  with  ice  con¬ 
centration  in  both  drift  phases  suggests  two  possibilities. 
First,  noise  is  parametrically  dependent  upon  ice  concentra¬ 
tion.  That  is,  as  ice  concentration  increases  so  does  the  den¬ 
sity  of  potential  surface  sound  sources.  Some  external  forc¬ 
ing,  perhaps  surface  gravity  wave  interaction  with  sea  ice, 
may  then  generate  low-frequency  sound.  Second,  the  actual 
convergence  of  the  ice  field  may  cause  floes  to  collide  and 
produce  low-frequency  sound.  For  example,  in  addition  to 
wave-induced  floe  collision,  tapping  vibration  is  also  plausi¬ 
ble  for  horizontal  collisions  resulting  from  differential  wind 
and  current  forcing.  The  differential  velocity  between  floe 
centers  is  roughly  s£div ,  where  s  =  DC  ~  1/2  is  the  mean  spac- 


TABLE  V.  Cross  correlation  between  ambient  noise  (rms  pressure)  in  the 
25-  to  50-Hz  band  and  ice  field  descriptors.  Environmental  time  series  are 
held  fixed.  (Data  sampled  at  6-min  intervals.)  [Peak  correlation  coeffi¬ 
cient,  time  lag  (h) .  ] 


1st  drift 

\P\ 

2nd  drift 

1^1  • 

^div 

-0.25,  -2.0 

^12 

0.25,  1.4 

^XX  ^yy 

0.24,  0 

c 

0.43,  -  1.0 

0.77,  -2.6 

^ xx  ^yy 

-0.15,  -  1.6 

"abs 

0.29,  0 

l^abs  I 

-0.15,  -1.0 

(e„iv  )2  +  (fii  )2 

-0.11,  0 
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ing  between  floe  centers  and  D  =  100  m  is  a  typical  ice  floe 
size  for  the  MIZ.  Low-frequency  sound  would  be  radiated  at 
/=  Dediv/8jC.  For  the  given  parameters,  and  for  large  va¬ 
lues  of  C~0.7  and  ediv<1.7x  10  _4  s~ 1  as  measured  during 
compact  conditions, /is  less  than  or  roughly  equal  to  20  Hz 
and  is  possibly  in  the  range  of  interest. 

III.  HIGH  NOISE  EPISODES  OCCUR  DURING  WEAK 
SURFACE  WAVE  FORCING 

During  off-ice  winds,  surface  wave  activity  was  weak 
and  episodes  of  high  noise  were  sometimes  observed.  These 
episodes  were  of  roughly  3-  to  6-h  duration,  1-  to  2-h  e-fold¬ 
ing  time,  recurred  at  12-h  periods,  but  were  present  only  7% 
of  the  time.  Portions  of  the  noise  record  specifically  related 
to  these  episodes  appear  in  Fig  7. 

We  cross  correlated  these  portions  of  the  noise  record 
with  environmental  forcing  functions,  as  summarized  in  Ta¬ 
ble  VI.  Although  an  argument  can  be  made  for  wind  stress, 
current  stress,  and  drift  acceleration,  the  best  correlate  of  the 
noise  is  consistently  drift  speed.  This  is  not  surprising  since 
drift  speed  records  display  substantial  12-h  periodicity  and 
have  episodes  corresponding  to  those  found  in  the  noise,  as 
shown  in  Fig.  7.  The  12-h  periodicity  in  drift  speed  is  likely 
not  caused  by  tidal  current  since  anomalously  strong  diur¬ 
nal,  rather  than  semidiurnal,  tides  characterize  the  MIZ 
near  the  site  of  the  experiment.26  Also,  current  measure¬ 
ments  at  50-m  depth,  where  momentum  transferred  from 
the  surface  is  minimal,  show  poor  correlation  with  the  noise, 
suggesting  tidal  currents  are  likely  unrelated  to  it.  The  semi¬ 
diurnal  periodicity  in  drift  speed  is  more  likely  related  to 
wind  driven  inertial  oscillations  of  the  ice  field.  Inertial  oscil¬ 
lations  can  cause  oscillating  convergence  and  divergence  of 
the  ice  field  which  could  result  in  high  noise  during  periods 
of  increased  concentration. 

Finer  spatial  and  temporal  sampling  was  available  for  a 
particular  noise  episode.  Over  a  roughly  1 2-h  period  at  about 


TABLE  VI.  Cross  correlation  between  noise  in  the  25-  to  50-Hz  band,  spe¬ 
cifically  related  to  episodes,  and  environmental  forcing  functions.  Environ¬ 
mental  time  series  are  held  fixed.  (Data  sampled  at  1/2  h  intervals.)  [Peak 
correlation  coefficient,  time  lag  (h).] 


1st  drift 

1*1 

2nd  drift 

1*1 

Kl 

0.30,  0 

0.70,  -0.50 

kJ 

0.60,  -0.5 

k«  +rj. 

0.18,  -0.5 

k-rj 

0.04,  -0.5 

|V,| 

0.72,  1.5 

0.77,  0.5 

\dVt/dt\ 

0.6,  -1.5 

\ap\ 

0.51,  -2.5 

Mr 

0.20,  -  1.5 

sw 

0,  0 

-0.18,  -1.6 

|VJ 

0.25,  2.5 

J171.5,  the  noise  records  of  seven  drifting  sensors,  averaged 
at  and  sampled  over  a  30-s  period,  are  shown  in  Fig.  8.  A 
chart  of  the  positions  of  each  sensor  during  the  episode  is 
shown  in  Fig.  9.  Comparing  these  positions  with  the  noise 
records,  we  observe  some  interesting  phenomena.  First,  sen¬ 
sors  furthest  to  the  southeast  measure  noise  of  greater  ampli¬ 
tude  during  the  episode,  suggesting  either  greater  density  of 
sources,  greater  source  strength,  or  greater  proximity  to  the 
sound  sources.  Second,  during  the  episode,  sensors  within 
about  1  km  of  each  other  measure  similar  temporal  varia¬ 
tions,  while  those  measured  by  sensors  at  separations  greater 
than  1  km  are  quite  different.  Finally,  the  episode  propagat¬ 
ed  through  the  array  of  sensors  at  roughly  2  to  5  km/h,  for  a 
mean  of  about  3  km/h  in  a  northwest  direction. 

This  propagation  speed  is  an  important  clue  to  the  epi¬ 
sode  mechanism.  It  propagates  much  too  slowly  to  be  an 
acoustic  wave.  It  is  also  too  slow  to  be  a  surface  gravity  wave 


FIG.  7.  Records  showing  high 
noise  episodes  during  weak 
surface  wave  forcing.  The  epi¬ 
sodes  are  of  3-  to  6-h  duration 
and  occur  during  off-ice  winds. 
Two  prominent  episodes  ap¬ 
pear  during  the  first  drift  phase 
(a)  and  are  separated  by  a  24-h 
period.  There  is  evidence  of  an 
episode  midway  between  them. 
Two  prominent  episodes  also 
appear  during  the  second  drift 
phase  (b)  and  are  separated  by 
12  h.  Corresponding  ice  drift 
speed  records  are  shown  in  (c) 
and  (d).  Their  significance  is 
discussed  in  the  text. 
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FIG.  8.  Records  of  rms  noise  in  the  25-  to  50-Hz  band  during  the  episode 
occurring  near  J171.5.  Corresponding  sensors  are  marked  to  the  right  of 
each  record.  Records  are  stacked  so  that  sensors  further  to  the  northwest 
are  lower.  The  episode  apparently  propagates  to  the  northwest  and  is  of 
smaller  amplitude  further  to  the  northwest. 


by  at  least  a  factor  of  5.  On  the  other  hand,  it  is  much  faster 
than  the  drift  speed  of  ice,  by  at  least  a  factor  of  3.  However, 
long  internal  waves  in  the  Fram  Strait  MIZ  can  travel  non- 
dispersively  with  speeds  in  the  range  0.90  <c<  6.8  km/h, 
according  to  Meunch  et  al.21  where 

and  where  g  is  the  acceleration  of  gravity,  H  is  the  mixed  or 
upper  layer  thickness  defined  by  the  pycnocline,  p  is  the 
mean  density  of  sea  water  and  Ap  is  the  density  jump 
between  upper  and  lower  ocean  layers.  ( Long  internal  waves 
have  wave  numbers  k{  such  that  The  noise  epi¬ 

sode’s  propagation  speed  is  within  the  range  of  long  internal 
waves.  Also,  the  speeds  are  consistent  with  density  jump  and 
mixed  layer  thickness  measured  during  MIZEX  84. 28  Ther¬ 
efore,  we  conclude  that  this  particular  noise  episode  is  affect¬ 
ed  by  the  passage  of  long  internal  waves. 

Meunch  et  al.  speculate  that  long  internal  waves  may  be 
generated  at  the  ice  edge  during  off-ice  wind  forcing,  and  this 
coincides  with  the  conditions  during  which  we  observed 
noise  episodes.  They  further  speculate  that  these  internal 
waves  may  concentrate  the  ice  field  into  bands  of  high  con¬ 
centration  due  to  alternating  convergence  and  divergence 
zones  created  at  the  sea  surface  along  the  wave’s  path.29 

We  speculate  that  the  propagating  noise  episode  is  due 
to  the  advection  of  bands  of  highly  concentrated  ice  by  a  long 
internal  wave.  We  further  speculate  that  the  other  noise  epi¬ 
sodes  observed  were  also  related  to  the  formation  of  bands  of 
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FIG.  9.  Chart  of  the  propagation  of  the  J171.5  noise  episode  through  the 
acoustic  array.  The  relative  arrival  time  of  the  episode  was  computed  for 
each  sensor  by  cross  correlating  each  noise  record,  shown  in  Fig.  8,  with 
that  of  the  first  sensor  to  detect  the  episode,  D19.  The  triplets  (sensor  name, 
correlation  coefficient,  time  lag  in  hours)  are  centered  at  the  position  of 
each  sensor  during  peak  correlation.  Sensors  deviated  from  the  relative  po¬ 
sitions  shown  by  less  than  1  km  during  the  total  episodic  time  of  roughly  3  h. 
Contours  have  been  drawn  to  show  positions  of  equal  time  lag  of  the  epi¬ 
sode.  The  episode  traveled  at  roughly  3  km/h  to  the  northwest. 


highly  concentrated  ice  by  internal  waves.  Additional  sup¬ 
port  for  this  appears  in  Fig.  10.  Bands  of  high  ice  concentra¬ 
tion  appear  in  the  vicinity  of  the  acoustic  array  at  the  same 
time  of  day  as  the  noise  episodes  which  occurred  two  and 
three  days  earlier,  and,  as  Fig.  4  shows,  during  which  the 
surface  waves  remain  equally  weak.  The  spacing,  width  and 
orientation  of  the  bands  suggest  the  propagation  of  a  long 
internal  wave  in  roughly  the  same  direction  as  the  noise  epi¬ 
sode. 

The  weaker  correlation  between  ice  concentration  and 
noise  found  during  the  first  drift  phase,  in  comparison  with 
the  second,  is  likely  due  to  the  inability  of  the  strain  array  to 
resolve  the  smaller  scale  ice  bands  which  we  speculate  to 
have  affected  the  noise  episodes. 

IV.  SUMMARY  AND  DISCUSSION 

We  have  studied  low-frequency  ambient  noise  records 
measured  under  drifting  ice  floes  of  the  Fram  Strait  MIZ 
under  early  summer  conditions.  We  find  that  long  term  var¬ 
iations  of  the  noise  are  strongly  related  to  the  amplitude  of 
open  ocean  surface  gravity  waves  incident  on  the  ice  edge.  In 
this,  we  confirm  the  conclusion  of  earlier  researchers,2  who 
attributed  MIZ  noise  to  surface  gravity  waves  interacting 
with  the  floes.  We  also  find  that  these  long  term  variations  of 
the  noise  are  related  to  ice  concentration  which,  in  turn,  is 
determined  primarily  by  wind-induced  free-drift  of  the  ice. 
While  ice  drift  exhibits  semidiurnal  inertial  oscillations, 
noise  associated  with  surface  wave  forcing  has  at  most  weak 
semidiurnal  components.  The  noise  under  surface  wave 
forcing  is,  instead,  created  and  integrated  over  large  spatial 
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FIG.  10.  A  side-looking  airborne  radar  image  at 
19  GHz  taken  during  a  NOAA  P3  flight  over 
the  ice-edge  at  about  J 1 74. 5. 25  The  swath  width 
is  roughly  1 1  km  and  the  distance  between  the 
arrows  on  the  horizontal  is  19  km.  The  Kvitb- 
jom  is  located  at  80.26°  N  7.4°  E,  which  is  just 
off  the  image  to  the  east.  The  white  specks  are 
ice  floes  of  roughly  100  to  200  m  in  size,  and  the 
ice  concentration  is  between  0.2  to  0.4  except 
where  bands  of  roughly  1-  to  2-km  spacing,  1- 
km  width  and  2-  to  3-km  length  appear  centered 
at  80.3°  N  between  6.2°  and  6.5°  E,  where  much 
higher  concentrations  are  found.  (The  vertical 
line  at  about  7.3°  E  is  due  to  splicing  of  the  im¬ 
age.  It  does  not  represent  an  actual  ice  band. ) 


(>5  km)  and  temporal  (>3  h)  scales,  where  the  numerical 
values  come  from  measured  correlations. 

Shorter  spatial  and  temporal  scales  can  also  be  observed 
in  the  data,  and  while  the  noise  at  such  scales  occurs  infre¬ 
quently  and  is  therefore  of  less  practical  importance,  its 
study  leads  to  confirmation  of  important  physical  ideas. 
Brief  episodes  of  high  noise  were  detected  during  off-ice 
winds  and  correspondingly  weak  surface  wave  amplitudes. 
These  episodes  are  of  roughly  1  h  e-folding  time,  of  semidiur¬ 
nal  periodicity,  and  are  strongly  related  to  ice  drift  speed, 
which  also  is  semidiurnal.  Because  such  episodes  propagate 
at  a  speed  consistent  with  internal  waves,  we  speculate  that 
they  are  additionally  affected  by  the  formation  and  advec- 
tion  of  bands  of  high  ice  concentration  by  internal  waves 
formed  at  the  ice  edge  during  off-ice  winds.  Thus  ice  concen¬ 
tration  has  a  twofold  effect  in  episodic  behavior,  first  via 
semidiurnal  inertial  oscillations  of  the  ice  floes  in  free  drift, 
and  second  via  ice  banding  driven  by  long  internal  waves. 
Such  ice  concentration  effects  combine  to  give  relatively 
shorter  spatial  and  temporal  scales  to  the  noise  than  that 
caused  by  on-ice  surface  gravity  waves.  When  the  latter  are 
present,  the  wind  is  expected  to  be  on  ice,  and  the  presumed 
generation  mechanism  for  long  internal  waves  would  be  ab¬ 
sent,  as  would  the  episodes  in  the  form  observed. 

We  conclude  for  both  long  and  short  scales  that  ice  con¬ 
centration  is  an  important  parameter  of  MIZ  noise.  It  deter¬ 
mines  not  only  the  spatial  and  temporal  scales,  but  also  indi¬ 
rectly  the  magnitude  of  the  noise.  In  this  sense,  earlier 
observations  of  MIZ  hotspots  attributed  to  ice-edge  ed¬ 
dies3,4  are  likely  of  the  same  class  as  the  observations  here, 
the  eddy  merely  imposing  spatial  and  temporal  scales  char¬ 
acteristic  of  its  dynamics,  in  contrast  to  the  ice-edge  dyna¬ 
mics  observed  in  the  present  study. 

It  should  be  no  surprise  that  a  stress  parameter,  such  as 
imposed  by  incident  surface  gravity  waves,  and  a  strain  para¬ 
meter,  such  as  ice  concentration,  jointly  relate  to  MIZ  am¬ 
bient  noise.  In  the  MIZ,  stress  and  strain  are  uncoupled;  no 
constitutive  law  for  the  field  can  exist  under  free-drift  condi¬ 
tions.  Further,  we  have  focused  on  the  one  best  correlate  for 
each  stress  or  strain  parameter.  It  is  likely  that  multivariate 
analyses  in  each  would  provide  firmer  knowledge  of  MIZ 
ambient  noise,  but  some  useful  physical  knowledge  can  non¬ 
etheless  be  extracted  from  the  narrower  focus  in  our  study. 

From  data  in  our  study,  we  can  assess  the  plausibility,  at 
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least  to  first  order,  of  the  actual  mechanisms  by  which  low- 
frequency  sound  is  generated  in  the  MIZ.  Because  surface 
gravity  wave  amplitude  is  the  best  noise  correlate  for  large 
spatial/temporal  scales,  physical  mechanisms  should  be  in¬ 
voked  involving  it.  First,  incoming  waves  can  cause  flexural 
fracture  in  ice  floes,  and  second,  they  can  cause  relative  verti¬ 
cal  and  horizontal  motion  between  adjacent  floes.  The  first 
radiates  noise  directly  from  each  fracture,  and  subsequently 
from  unloading  motion  as  each  fractured  floe  returns  to  its 
undeformed  state.20,22  The  second  radiates  noise  as  asperi¬ 
ties  on  edges  of  adjacent  floes  in  contact  break  or  vibrate  in 
response  to  relative  floe  motion.  For  the  present  study,  sur¬ 
face  wave  amplitudes  within  several  kilometers  of  the  ice- 
edge  were  sufficient  to  cause  floe  fracture  and  asperity  break¬ 
ing  or  vibration.  The  frequency  range  within  which  flexural 
fracture  and  unloading  motion  creates  noise  is  estimated  to 
easily  cover  the  low-frequency  range  of  interest.  Asperity 
vibration  can  cover  portions  of  this  range  at  the  lowest  end, 
but  asperity  breaking  is  too  high  in  frequency.  So  we  nomin¬ 
ate  flexural  fracture  and  asperity  vibration  as  the  plausible 
mechanisms.  Flexural  fracture,  however,  is  much  more  like¬ 
ly,  not  only  because  of  its  broader  frequency  coverage,  but 
also  because  much  larger  energy  is  invested  in  a  floe  to  frac¬ 
ture  rather  than  to  vibrate  it,  with  noise  radiation  therefore 
larger. 

Small  scale  variations  of  the  noise,  as  in  episodes,  are 
unrelated  to  surface  waves.  It  is  possible  that  floe  collision  in 
ice  convergence  is  responsible,  but  we  estimate  that  this  me¬ 
chanism  covers  the  range  only  for/ <  20  Hz  (see  comment  at 
the  end  of  Sec.  II B).  The  restricted  frequency  range,  plus 
the  low  overall  correlation  of  ediv  with  the  noise  (see  Table 
V),  prevents  us  from  labeling  this  mechanism  as  plausible. 
While  the  source  mechanism  thus  remains  unsettled  for 
small  scale  variations,  we  do  not  have  a  better  candidate  than 
floe  collision  in  ice  convergence. 
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APPENDIX 

To  obtain  quantitative  measurement  of  temporal  scales 
and  the  linear  relationship  between  time  series,  a  cross-cor¬ 
relation  algorithm  was  devised  that  requires  no  interpolation 
of  data  gaps  and  so  introduces  no  hypothetical  information. 
Time  series  had  only  to  be  resampled  to  a  consistent  sample 
period  upon  which  gaps  were  reintroduced.  For  each  lag, 
discrete  correlation  coefficients  could  be  computed  by  sum¬ 
ming  only  overlapping  samples  and  so  ignoring  gap-gap  and 
gap-sample  overlaps  in  the  convolution.  Time  series  in  each 
lag  window  were  locally  “zero-meaned”  and  normalized  by 
the  square  root  of  the  variance  product  of  each  time  series 
using  only  nonzero  overlap  samples  in  the  local  lag  window. 
Even  a  time  series  with  no  gaps  would  require  a  moving 
normalization  window  to  account  for  the  nonstationarity  in¬ 
herent  in  data  streams  of  finite  length.  The  scheme  is  symbo¬ 
lically  described  below. 

The  ideal  time  series  with  no  data  gaps  are  x[n]  and 
y[n],  of  lengths  Nx  and  Ny,  respectively.  Time  series  8x[n] 
and  8y[n]  are  zero  for  data  gaps  at  n  such  that  8X  [n]x[n] 
and  8y  [n]y[n\  are  the  actual  time  series  with  data  gaps. 
The  cross-correlation  function  is  then 

Rxy(m) 

^  ^  (*[n]  - x[m])(y[a ]  -y[m])8x[n]8y[a] 

«  =  o  ax[m]ay[m]N  [m] 

for 

0<m<Nx  +  Ny  —  1,  a  =  n  —  [m  —  (Ny  —  1)], 
where 

x[m]=N~'[m]  J  x[n]8x[n]8y[a], 

n  =  0 

y[m]  =N~1[m]  £  y[a\8x[n\8y[a\, 

n  =  0 

oiim]  =N -l[m]  ^  x2 3[n]8x[n]8y[a]  -  (x[m])2, 

n  =  0 

ojlm]  =N  -l[m]  ^  f[a]8x[n]8y[a]  -  ( y[m ])2, 

n  =  0 

N[m]=  ^  8x[n]8y[a\. 

n  =  0 

The  same  procedure  was  used  for  all  correlations  in  this 
paper. 
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